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We analyze the exact ground state of XXZ zigzag spin chain with apphed magnetic field and find 
the quantum critical surface. Using the theorem of positive semi-definite matrix, we can prove that 
the ground states for a specific region, are fully polarized state and one magnon states. With Bethe 
ansatz, we argue that this is the quantum critical surface in all cases. A first in the literature, we 
derive the analytical expression of quantum critical surface for superconducting quantum dots array 
in the presence of gate voltage. 
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Low dimensional spin systems has been an active field of 
research [l| for many year. Bethe ansatz has been used 
widely to solve many new families of integrable Hamil- 
tonians of one dimensional systems 043- A quantum 
inverse scattering method based on Bethe ansatz has 
also been developed for dealing with the two-dimensional 
models f^. Another method which also enables one to 
find the exact solutions of various models is the matrix 
product method 0413 • For example, it is related to 
the AflJcck-Kennedy-Lieb-Tasaki(AKLT) state ^E^. In 
this work we are going to apply these two methods to the 
zigzag spin chain model. 

A zigzag spin chain is an one dimensional system of spins 
with exchange interactions with their nearest (NN) and 
next-nearest neighbors (NNN) , denoted by Ji and J2 re- 
spectively. It is a frustrated system if J1J2 < or both 
of them are positive. In a zigzag spin chain the situa- 
tion is not clear. Majumdar and Ghosh 12] found that 
the degenerate dimer states in which neighboring sites 
form singlets (MG state) are the exa ct g round states at 
^ = 2. At ^ = -4, Hanada et.al. ^ found that the 
uniformly distributed resonate valence bond (UDRVB) 
state and the fully polarized (FP) state are degenerate 
ground states. There are also a few other important the- 
oretical works in zigzag spin chain [14,-J/7j. The zigzag 
spin chain can be realized in physical systems, such as 
the oxides, LiaCuOa (multiferroic) ^ and SrCuOa 
In this letter, we find a quantum critical surface on which 
the FP state is degenerate with the one-magnon states, 
along with an application to superconducting quantum 
dots system. Our calculation is based on the concept 
that the local ground states in a matrix product form 
can be extended to global ground states. 
Model Hamiltonian: We begin with the Hamiltonian 
of a zigzag spin chain with applied magnetic field. 
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where Ji and J2 are the NN and, NNN interaction, 
Ai(A2) is the anisotropy in z axes for NN(NNN) inter- 
action and i is the label of the lattice site. Furthermore, 
we use the periodic boundary condition and hence we 
have the coupling between sjv-i,sjv and si,S2. We can 
dissect the Hamiltonian into many local Hamiltonians. 
Each local Hamiltonian contains three spins. Thus, the 
original Hamiltonian can be written as 
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-|-sf_,_2). i+2 is the local Hamiltonian 

of exchange interaction, and h^i^i is the that of ap- 
plied magnetic field. We have set 6 = -^,f = ^ and 
made scaling J2 = 1- In we have introduced a free 

parameter x to facilitate later calculation. Note that the 
Hamiltonian is independent of the value of x if we apply 
the periodic boundary condition. 

The quantum critical surface from local site 
Hamiltonian: The local Hamiltonian /ii23, which serves 
as our starting point, can be diagonalize and the eigen- 
values in terms of parameters x, Ai^2 and / are 
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where 



in a more compact form 
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We found that Eq{Ei) corresponds to the state 
(lttt))i E2 and E^± correspond to one spin-up 
states and and £'4-1- to the one spin-down states. 
In order to see quantum critical points for entire system, 
we seek in the spectrum the level crossing between FP 
state I 4,4,4,} and other eigen state. We need at least two 
other states because for an infinite spin chain with space 
inversion symmetry any other states must be two-fold 
degenerate. Inspecting the eigenvalues, we found that 
Eq,E2 and £'5- are three lowest possible energies. By 
requiring Eq — E2 = £5-, we get the following relations 
^ = ^JiitttfV & = 3(2 + 2A2 + Ai/))(2 -f x) or for 
indication of the quantum critical surface 
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The correspond energy of these three degenerate state is 
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Matrix formulation of the global eigen state: Here, 
we generalize the results of local Hamiltonian for arbi- 
trary number of spins. The global Hamiltonian of N 
spins can be obtained by the direct sum of local Hamil- 
tonian. 
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where each basic element h is just the local Hamiltonian 
hi.i+i,i+2 and 1m is the identity matrix of rank AI . The 
eigen states of an infinite spin chain can be also found by 
adding more spins to eigen states of the local Hamilto- 
nian. Consider the matrix 



(4) 



representing the spin state at site j. The three degen- 
erate ground states of the local Hamiltonian hjj+ij+2 
with energies Eq, E2 and E^^ can be expressed as the 
matrix mjmj^imj^2- In fact, the matrix elements are 
the linear combinations of these three eigen states. It 
is clear that Af„ — Hj^^TOj gives the states of n spins. 
hjj+ij+2, operates only on the product mjmj+imj+2 
(whose elements are its eigen states,) without affecting 
other matrices. Hence, the operation of Hamiltonian will 
not alter the form of M„ and its matrix elements are lin- 
ear combinations of the eigen states. We can write Af„ 
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where |0 = |0),|<^:;) = E;^! lj) and = 
E;=i/«-jIj), with |0) = I UU = and 

the coefhcients fj satisfying the relations /o = 0, /i = 
1' fj + 2 /j+i + /j+2 ~ Note that fj — cos{jk) where 
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satisfies above equations. The first eigen state is the 
fully polarized (FP) state and the other two can be ex- 
pressed as the linear combination of two independent 



one-magnon states: 



) - -^T"- e 



ikj + 



s,'|0), 10=) 



^ Ej'=i where j is the site position. The 

periodic boundary condition implies that Nk = 2rmT. 
Thus we have constructed the eigen states for a zigzag 
spin chain of arbitrary length at the surface Eq.([2]). 
Proof of the quantum critical surface by positive 
semi-definite matrix theorem: We use the theorem of 
positive semi-definite matrix to show that for spin chain 
of arbitrary length, the above mentioned three degen- 
erate eigen states are actually the global ground state 
in certain region. The theorem of positive semi-definite 
matrix is the following: The necessary and sufficient con- 
dition for a real symmetric matrix A to be positive semi- 
definite is x'^Ax > for all real vectors x. If M and N 
are positive semi-definite, then the sum M -\- N, and the 
direct sum M (B N, and direct product M <Si N are also 
positive semi-definite. 

The FP state and the one-magnon states are the de- 
generate ground states, if the matrix (/ijj+ij+2 — EqI), 
with Eq being their energy, is a positive semi-definite 
matrix. We begin the proof with a local Hamiltonian 
with three spins. By setting x — 2(2'+2A2+a\/) ' Ecl-© 
and writing the Hamiltonian hi^i^i^i^2 — Eq in a real and 
symmetric matrix form, we proceed to prove it is a posi- 
tive semi-definite matrix. The value {hi^i+i^i+2 — Eo)x 
where x = {xi,X2,X3,X4,X5,xq,xj,xs)'^ is any real col- 
umn matrix, is 
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{hj,j+i,j+2 — EqI) is a positive semi-definite matrix pro- 
vided two inequalities below are satisfied simultaneously 
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Since all the local Hamiltonian {hjj+ij+2—Eoi) are pos- 
itive semi-definite matrices, so is /1123 (8 12 + 12 "Xi ^234 and 
H — Eq. Thus we conclude that in the region given by in- 
equalities ([7]) , the degenerate ground state is the FP state 
and one magnon state with k in Eq.®. The degenerate 
ground states are I^n) = |0) , \^%) - E™=i e^""s+ |0) 
I'^Af) — J2m=i ^~''^"^^m We will give argument below 
that even if outside of the region given by inequalities 
(O, they are still the degenerate ground states. In the 
mean time we show the quantum critical surfaces given 
by Eq.® for the case Ai = A2 = d in Figdja) and the 
case Ai = d and A2 = in Figlljb). In the region above 
the surfaces only the FP state is the ground state. There 
are several interesting features. Firstly, larger Ai and 
A2, requires larger magnetic field b. This means that the 
system is made difficult to be polarized by anisotropy. 
On the other hand, negative /, ferromagnetic coupling, 
reduces the value of h as expected. Eq.® is valid for 
/ > —4. Beyond which only the FP state is the ground 
state. Quantitatively, we see a linear dependence of b 
on Ai and A2. This indicates that treating the coupling 
SiS^ with the mean-field approximation is accurate at 
least in the region close to that given by Eq.®. 

Analytical solutions outside the region of Eq.([7|)): 

Beyond the regions given by Q , the global Hamiltonian 
is still likely to be positive semi-definite. We apply Bethe 
ansatz to analyze this region. Our approach is the fol- 
lowing. We write down the Bethe ansatz solution for n 
up-spins state, and the constraints of the coefficients. For 
a given h and a given /, we find the state with lowest en- 
ergy. We are able to show that for n < y, the energy of 
above mentioned state is greater or equal to the energy 
of FP state on the surface given by Eq.®. We begin by 
writing the state as the linear combinations of n spin-up 
states: 

IV') = o.xux^,...,x^st^st^...stjO) (8) 
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where axi,x2,...,xn is symmetric under exchanges of spin 
indices. Their constraints can be obtained by requiring 

to be the eigen state of the Hamiltonian. 
In this case there are many constraints. But one of them 
has the property Xi+i > Xi + 2 with all possible i. This 
means that all the up-spins are separated by at least two 
down-spins, or no two up-spins are NN or NNN. This 
relation is 
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FIG. 1; The quantum critical surface in the case (a) Ai = 
A2 = d, and (b) Ai = d, A2 = 0. 

With the same method for two-magnon states, we let 
(ixi,X2,...,X„ = J2p cxp[i X]Li kpjXj]^ where P means 
permutation of spin sites j. The equation becomes 

n n 

2{E-Eq) = 2n(6-A2-Ai/)+2/^cosfc,+2^cos2fc,. 

i=l i=l 

(9) 

If ki's were independent variables, then the lowest en- 
ergy would have been obtained by varying them inde- 
pendently. The result is 

ki ^ k ^ ±cos-\-^). (10) 

For E ~ Eq and Eq. pH]) . we again reach the relation in 
Eq. ®. Since fc^'s have other constraints, the energy we 
thus get is the lower bound. Hence, we conclude that for 
1 < < y J the FP state and one-magnon states have the 
lowest energy with the magnetic field given by Eq. ®. 
It is plausible to assume that this is true for arbitrary 
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value of n. Bethe ansatz gives the excitation energy of 
magnons as AE — E — Eo = l + ^+ f cos k + cos 2k. 
Interestingly, it does not depend on the anisotropy Ai 
and A2. It vanishes when Eg.© is satisfied. Physically, 
this means that the magnons soften at certain wave vec- 
tors, depending on the value of /. This phenomena can 
be detected by neutron scattering or susceptibility mea- 
surement. 

Application to superconducting quantum dot sys- 
tem: Here we derive the analytical expression of the 
quantum critical surface of a superconducting quantum 
dots (SQD) array with NN and NNN tunneling and 
Coulomb interaction in the presence of gate voltage in 
every superconducting dot. The quantum critical sur- 
face separate the charge density wave phase of lower 
commensurability {n = 1/2), which corresponds to a 
anti-ferromagnetic phase of a spin chain; from a higher 
commensurability state (n = 1, or integer numbers of 
Cooper pair in every dot,) which correspondence to the 
FP phase. Therefore, we use the relation between b 
and / evaluated in the previous section to study the 
quantum critical surface between these two phases. The 
quantum fluctuations are controlled by the system pa- 
rameters like charging energy of SQD and the Joseph- 
son coupling (Ej). In the SQD model [21|, we con- 
sider the finite range Ej and Coulomb interaction which 
have experimental support 22-2^. The model Hamil- 
tonian of our study consists of different interactions, 
H = Hji + Hj2 + Heco + Heci + Hec2- We recast 
our basic Hamiltonians in the spin language around the 
charge degeneracy point. Hji = 2EjiJ2i{Si'^Si+i^ + 

h.c), Hj2 = 2E.J2Y.^{S^^S,+2- + h.c), HeCO = 

^J:^{2S,'' + h)^ « 2hEcoE^^S,'' , Heci = 
4-Bzi J:^ S,+i^, Hec2 = 4Ez2 J:^ S,+2^ ■ whcrc 
Hji and Hj2 arc Josephson energy Hamiltonians respec- 
tively for NN and NNN Josephson tunneling between the 
SQD and Heco, Heci, and Hec2 are respectively the 
Hamiltonians for on-site, NN and NNN charging ener- 
gies of the SQD. h = ^"+1-^ is a parameter which al- 
lows tuning the system to a degeneracy point by means 
of gate voltage. When Eco » Ej, sequential tun- 
nelling of Cooper pair across the SQD is not a ener- 
getically favorable process. One must consider the co- 
tunneling effect, i.e., the higher order expansion in 
which reduce the intersite Coulomb charging energy and 

M |27|. 



also increase the NNN Ej [21|, |26|, l27| . This results in 
new tunneling energies: Ej2 Ej2 + Eji^/2Ecq, and 
Ezi Ezi — iEji^ /\&Eco- 

Now we present the analytical expression of quantum 
critical surface for the presence and absence of co- 
tunneling effect, following the relation h and / in Eq.®. 
The analytical expression of quantum critical surface in 
the absence of co-tunneling effect are the following 
2hEca = 2Sj2 + 4(^,1 + E,2) + ""''^ 



2hEco = 2Ej2 + 4{E,2 + E,i) + ^ + ^ + ^.We 
have already calculated the quantum critical surface on 
which FP state and one magnon state are degenerate. 
Here the one vacant site of Cooper pair in the SQD lat- 
tice array corresponds to one magnon state and Mott 
insulating state with integer numbers of Cooper pairs in 
each lattice site corresponds to the FP state. The ap- 
plied gate voltage on each SQD of the array corresponds 
to the applied magnetic field. The charging energy cor- 
responds to the anisotropy of the exchange interaction of 
the spin system. For large gate voltage, the SQD is in 
the insulating state. Interestingly, the charging energy 
tends to reverse this trend by destabilizing the state of 
integer numbers of Cooper pairs in each lattice site and 
makes the system conducting. The analytical expression 
for the quantum critical surface for the SQD array is the 
first in the literature. 

Conclusions: We have found the quantum critical sur- 
face of a zigzag spin chain under magnetic field by using 
the theorem of positive semi-definite matrix and Bethe- 
ansatz methods. We have also find the quantum critical 
surface of superconducting quantum dots system by us- 
ing our analytical methods. 
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while that in the presence of co-tunneling effect is 
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